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Abstract — The optimal bit-wise demodulator for Af-ary pulse 
amplitude modulation (PAM) over the additive white Gaussian 
noise channel is analyzed in terms of uncoded bit-error rate 
(BER). New closed-form BER expressions for 4-PAM with any 
labeling are developed. Moreover, closed-form BER expressions 
for 11 out of 23 possible bit patterns for 8-PAM are presented, 
which enable us to obtain the BER for 8-PAM with some of 
the most popular labelings, including the binary reflected Gray 
code and the natural binary code. Numerical results show that, 
regardless of the labeling, there is no difference between the 
optimal demodulator and the symbol-wise demodulator for any 
BER of practical interest (below 0.1). 

Index Terms — Additive white Gaussian noise channel, binary 
reflected Gray code, bit error probability, bit-interleaved coded 
modulation, demapper, demodulator, LLRs, logarithmic likeli- 
hood ratio, pulse-amplitude modulation, uncoded transmission. 



I. Introduction and Motivation 

Current wireless communication systems are based on the 
bit-interleaved coded modulation (BICM) paradigm introduced 
in |1| and later studied in |2|, |3|. One key element in 
these systems is the demodulator which calculates logarith- 
mic likelihood ratios (LLR, also known as L-values) for the 
received bits, which are then passed to the channel decoder 
The calculation of L-values is crucial in many other coded 
systems. The coded performance analysis of BICM systems is 
generally not straightforward and is usually carried out either 
numerically by Monte-Carlo simulation or in terms of lower 
and upper bounds [2, Sec. 4], [3, Ch. 4]. In this paper, we 
analyze the uncoded performance of bit-wise demodulators 
over the additive white Gaussian noise (AWGN) channel. 

The optimal bit-wise demodulator (BD) minimizing the 
BER implies the calculation of (exact) L-values for the re- 
ceived bits. The uncoded performance of such a demodulator 
has been studied in Q, where closed-form expressions for the 
BER for 4-PAM with the binary reflected Gray code (BRGC) 
im-Q are presented. Due to the complexity of the BD, the 
calculation of L-values in practical systems is usually done 
based on the so-called max-log approximation fS] eq. (5)], 
in eq. (1)]. We call this demodulator the approximate BD 
(ABD). The ABD is equivalent to the symbol detector in terms 
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of uncoded BER |10, Sec. IV-A], whose performance is well 
documented in literature, e.g., ill. Ch. 5], lil2, Ch. 10], [61, 
ifTsl - lfTsl and references therein. 

It is well known that the uncoded BER of one-dimensional 
constellation can be expressed as a sum of Gaussian Q- 
functions, cf. ifTTl Ch. 5], lfT2l Ch. 10] and references therein. 
The arguments of the Q-functions depend on the points that 
separate the decision regions associated to different bits. We 
refer to these points as thresholds. The computation of the 
thresholds for the BD — the optimal bit-wise demodulator — is 
in general complicated and unknown. In this paper, however, 
we show that this problem can be solved analytically for 4- 
PAM and any labeling, extending the results presented in [4]. 
Moreover, we also analytically calculate the thresholds for 
8-PAM with some relevant labelings, including the BRGC, 
the natural binary labeling (NBC) |[T9] Sec. II-B], the folded 
binary code (FBC) [16] USl Sec. II-B], the binary semi- 
Gray code (BSGC) HH Sec. II-B], and the anti-Gray code 
(AGC) f20l. Numerical results show that optimal and subop- 
timal demodulators are different in terms of the BER only 
at a very low SNR. At BER below 0.1 there is no notable 
difference between them. 

The rest of the paper is organized as follows. In Sec. Ullwe 
introduce the notation convention, the system model, and the 
two demodulators. In Sec. |III] the BER analysis is presented. 
The patterns that form a labeling are studied in Sec. |IV] The 
threshold computation for the BD is shown in Sec. |V] and the 
numerical results in Sec. |VI] The conclusions are drawn in 
SeclVn] 



II. Preliminaries 

A. Notation Convention 

In this paper the following notation is used. Lowercase 
letters x denote real or complex scalars and boldface letters 
X denote a row vector of scalars. The complex conjugate 
of X is denoted by x* . Blackboard bold letters X denote 
matrices with elements Xi,j in the zth row and the jth column 
and (•)^ denotes transposition. Calligraphic capital letters X 
denote sets, where the set of real numbers is denoted by 
TZ. The binary complement of x € {0, 1} is denoted by 
X — 1 — X and its bipolar representation by x — 2x — 1. 
Binary addition (exclusive-OR) of two bits a and b is denoted 
by a 6. Random variables are denoted by capital letters X 
and probabilities by Pr{ }. The Gaussian Q-function is defined 
as Q{x) ^ (l/Vt^) XT exp(-iV2) dt. 
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B. System Model 

In this paper we analyze a system where a vector of binary 
data b = [61 , . . . , bm] is fed to a modulator. The modulator 
carries out a one-to-one mapping from b to one of the M 
constellation points x E X = {si, . . . , sm} for transmission 
over the physical channel, where M — 2'". We assume that 
si < S2 < ■ ■ ■ < Sm- The modulator is defined as the function 
$ : {0, 1}™ X. 

The modulator is defined by the constellation and its binary 
labeling. A binary labeling is specified by the matrix C = 
[cj, . . . , c|j]^ of dimensions M by m, where the ith row 
Ci = [ci^i, . ■ . ,Ci.m] is the binary label of the constellation 
point Si, i.e., <i>(ci) — Si. 

For PAM constellations, Si = —d{M — 2i + = 
1, . . . , M, where d = ^Ji/^Iv'P — 1) to normalize the constel- 
lation to unit average energy, i.e., Es = (1/Af) X)f=i *f — 1- 
We assume the bits to be independent and identically dis- 
tributed (i.i.d.) with Pr{i?j = u} = 0.5, Vj and u e {0,1}, 
and thus, the symbols are equiprobable, i.e., PrjJf — si} = 
1/M, Vi. 

In this paper we consider a discrete time memoryless 
AWGN channel with output y ^ x + ri, where x E X and the 
noise sample 77 is a zero-mean Gaussian random variable with 
variance No/2. The conditional probability density function 
(PDF) of the channel output given channel input is 



PY\x{y\x) 



(1) 



where the average signal to noise ratio (SNR) is defined as 
7 ^ E,/No = 1/No. 

The observation y is used by the demodulator to decide on 
the received binary sequence, i.e., to produce b = [61, ... , bm]- 
In this paper we consider two demodulators to obtain b from 
y, which are described in the next section. 

C. Demodulators 

The BD calculates (a posteriori) L-values for the m bits 
based on the observation y, i.e., 

Fi{B, = l\Y = y} 



log 



Ft{B, = 0\Y = y} 

Ex 



-7(y-a:)2 ' 



(2) 



(3) 



for j = l,...,m and Xj^ — {si G X : Ci,j = w,Vi}. To 
pass from (|2]l to Q Bayes' rule was used together with the 
i.i.d. assumption of the bits and the conditional PDF in ([T]). 
The BD uses the L-values in ^ to make a decision on the 
received bit according to the rule 



1 ifl,{y)>0, 
otherwise. 



(4) 



The implementation of the BD in its exact form (|3]l is 
complicated, especially for large constellations, as it requires 
calculation of the logarithm of a sum of exponentials. To 
overcome this problem, approximations are usually used in 
practice. The most common approximation is the so-called 



max-log approximation (log^j e^' « max^ Xi) lU eq. (3.2)], 
fSl, eq. (9)], m eq. (5)], flU] eq. (8)], which used in Q gives 



min {y 



min (y 



xf 



(5) 



The use of the max-log approximation transforms the nonlin- 
ear relationship (|3) into a piece-wise linear relationship dU, 
as previously shown in e.g., [22, Fig. 3], f2T, eqs. (11)-(14)]. 

The ABD is defined as the demodulator that applies the 
same decision rule (|4]l based on L-values calculated by (|5]l. 
As mentioned in [ lOj Sec. IV-A] the ABD is equivalent to the 
symbol detector in terms of uncoded BER. 

III. BER FOR One-Dimensional Constellations 
The BER for a given labeling C can be expressed as 



^ rn 



(6) 



where using the law of total probability, the BER for the jth 
bit position Pj = Pr{Bj ^ bj\Bj = bj] can be written as 



P. 



1 



M 



M ^ 

1=1 



(7) 



The BER for the jth bit position Pj depends only on the 
subconstellations Xj,Q and A'j 1 (cf. ©-dD), i.e., Pj is a 
function of the jth column of C. 

We define a bit pattern (or simply pattern) as a length- M 
binary vector p — [pi,...,pA/] e {0,1}*^ with Hamming 
weight A//2. The labeling C can now be defined by m 
patterns, each corresponding to one column of C. We index 
the patterns as with w being the decimal representation 
of the vector p, i.e., w = X]f=i 2*^^'pi. For example, for 
M — 4, the pattern [0,1,0,1] is indexed as pg. The BER for 
the labeling C does not depend on the order of its columns, 
and thus, the BER for the labeling C is fully determined by a 
set of TO patterns (indices) >V = {wi, . . . , Wm]- Based on the 
previous discussion, from now on we concentrate our analysis 
only on patterns (and not on labelings). 

To analyze the BER of a pattern (PBER), the observation 
space Tl, is split into two disjoint decision regions, i.e., — 
{?/ e 7^ : 6 = 0} and 3^1 = {?/ e 7^ : S = 1} such that 
3^0 U 3^1 = 7^. Using the definition of 3^o and 3^i, the PBER 
for the pattern p = [pi , . . . , pm] can be rewritten as 

1 

By expressing P as in (O, it is clear that the PBER in (|7]) can 
be calculated using the decision region only, as opposed to 
alternative approaches where (|8} is expressed in terms of the 
PDF of the L-values |24l eq. (19)]. 



A. Decision Thresholds 

One key element in the BER analysis presented in this paper 
is the decision thresholds. Decision thresholds (or simply 
thresholds) for a given pattern p are the points that separate the 
decision regions for zeros and ones, and thus, they determine 
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Fig. 1. L-values in j3) vs. the received signal for different 7 for 8-PAM and the BD for pgg = [0, 1,0, 1, 0, 1, 0, 1]. Squares show the constellation points 
and crosses show the thresholds in Theorem |6] Gray and white areas indicate 3^0 ™d 3^1, resp. In (a) at 7 = 5.5 dB none of the seven thresholds is 
virtual. At 7 4.9 dB, shown in (b), the thresholds /3i and /32 as well as /Jg and fi-j become virtual, shown also in (c) for p = 4 dB. At 7 f» 2.2 dB, shown 
in (d), fio, and /Js merge with /34 and also become virtual. 



the PEER for the BD/ABD in For a given pattern p, 
we associate the threshold f3k E TZ to the bit pk when 
Pk 7^ Pk+1- Since there is no threshold j3k when pk — Pk+i, 
the number of thresholds is at most M — 1. The indices 
of the thresholds for the pattern p form a set of indices 
K,, with 1 < |/C| < M ~ 1. For example, the pattern 
P54 = [0, 0, 1, 1, 0, 1, 1, 0] has IC - {2, 4, 5, 7}. 

The thresholds for the ABD, which we denote by f3, are 
independent of 7 and placed at the midpoints between adjacent 
constellation points with different binary labels, which follows 
directly from (|5). On the other hand, the thresholds for the 
BD depend on 7. We denote these thresholds by /3, where to 
simplify the notation, the dependency on 7 is omitted. Fig. [T] 



shows the thresholds for the BD and their dependency on the 
SNR. This figure shows that some thresholds can merge at low 
SNR and seem to disappear To take this effect into account, 
we define virtual thresholds as follows. A threshold /3k is said 
to be virtual at 7 < 70 if it merges with another threshold f3k' 
at 7 = 7o (i.e., /3k = (3k' when 7 = 70) and does not exist at 
7 < 70. 

B. General Expression for One- dimensional Constellations 

The BER expression for the ABD and an il/-PAM con- 
stellation with any labeling is well known and can be found 
in l?! eq. (21)]. The PBER expression can easily be obtained 
in a similar way. The following theorem gives a general PBER 
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expression if none of the thresholds j3k is virtual. 

Theorem 1: The PEER of the BD or the ABD using an 
arbitrary one-dimensional constellation with a pattern p can 
be expressed as 

1 1 

i=l keIC 

where none of (3^ is virtual, and gi ^ G {±1} is 

9i,k = {Pk+i - Pfc)(l - 2p,;). (10) 

Proof: The proof is given in Appendix |A] □ 
The following theorem shows that Theorem [T] also holds 

when some of the thresholds become virtual, provided that 

their values are chosen properly. 

Theorem 2: If f3k is virtual for 7 < 70 because /3k = fik' 

at 7 = 70 for some k and k' > k, the PBER for the BD is 

given by ©-([lOll if /3fc = l3k', V7 < 70. 

Proof: We will show that the PBER in Q is not affected 

by the virtual threshold /3fe if (3k = f3k'- Let Si be the two 

terms in the inner sum in (|9]l associated to the thresholds Pk 

and (3k', i.e.. 

Si = 5i,/cQ((/3fc - 5^)1/27) + 5i,fc'Q((/3fc' - sOv^)- 

(11) 

Since (3k and f3k' are thresholds that merge, Pk+i ~ Pk' and 
Pk = Pk'+i must hold. Using these relations in ( fTOl i, we obtain 
5i,fc = —gi.k'i which used in (fTTT i gives 5^ = 0, Vi if (3k — 
(3k'. ^ □ 
Remark 1: Theorem |2] holds regardless of whether (3k' is 
virtual or not. If (3k' is not virtual for 7 < 70, (3k must be set 
to the value of (3k' for 7 < 70. If (3k' is virtual for 7 < 70, 
/?fc = (3k' can be set to any real value. 

C. BERfor the ABD and M-PAM 

For M-PAM and the ABD, Q can be expressed as a bit- 
wise version of |7, eq. (21)]: 

M-l 

^= M E «"Q((2"-1)o'v^), (12) 

where 

A/-1 

a„ = ^ (Pfc+i -Pk){^ ~ 2pk+i^n) 

k—n 

- {Pk+2-n -_Pfc+l-„)(l - 2pA;+l)- (13) 

One direct consequence of (fTSl i is that the vector a = 
[ai, . . . , flM-i] with a„ given by (fTsl i completely defines 
the performance of the ABD for 7\/-PAM and allows us to 
compare the performance of different patterns. From (ITZi . 
the PBER for high SNR is determined by the coefficient 
multiplying the Q-function with the smallest argument, i.e., 
ai. If two patterns have identical coefficients ai, the next 
coefficients 02 should be compared, and so on. 

Using ^ and the average BER for an Af-PAM with 



a labeling C can be expressed as 

P^-^H «„Q((2n - l)d^) , (14) 

n— 1 

where a. = [ax, ... , au-i] is the sum of vectors a for the m 
patterns used in C. The equation in ( fT4l i in fact corresponds 
to [7, eq. (21)], where the value of a„ is a scaled version of 
the so-called differential average distance spectrum 5{n, A), 
i.e., an = 2M5{n,\). 

IV. Bit Patterns 

For 4-PAM, the patterns pg = [0,1,0,1] and p^Q = 
[1,0,1,0] have identical PBER performance because of the 
symmetry of the constellation. To find all the patterns with 
different performance, we group patterns into classes, where 
all the patterns in the same class have identical PBER. 

In the following, we define two operations that can be 
applied to a pattern. A reflection of p is defined as p' = refl(p) 
with p'^ = pM+i-i for i — 1, . . . , M. For example, P27 = 
[0,0,0,1,1,0,1,1] = refi([l, 1,0, 1,1, 0,0,0]) = rc&{p^,e). 
An negation of p is defined as p' = nGg(p) with p'^ = pi 
for i ~ I, . . . , M. For example, = [0, 0, 1, 0, 0, 1, 1, 1] = 
neg([l, 1, 0, 1, 1, 0, 0, 0]) = neg(p2i6)- Both these functions 
are self-inverse, i.e., p = refl(refl(p)) and p = neg(nGg(p)), 
and they commute, i.e., refl(neg(p)) ~ neg(refl(p)). Note 
also that for any pattern p, we have that p 7^ neg(p). Using 
the previous definitions, we now define three special types of 
patterns that will be useful throughout this paper 

The pattern p is said to be reflected (RE) if rcfl(p) = p. For 
example, pgg — [0, 0, 1, 1, 1, 1, 0, 0] is an RE pattern. The pat- 
tern p is said to be anti-reflected (ARE) if ncg(rGfl(p)) = p. 
For example, P43 = [0, 0, 1, 0, 1, 0, 1, 1] is an ARE pattern. The 
pattern p is called asymmetric (ASY) if it is neither RE nor 
ARE. For example, P216 = [1,1,0,1,1,0,0,0] is an ASY 
pattern. 

From ©-(ISll, we note that the PBER is not affected by 
reflections and/or negation of the patterns, since the PBER is 
averaged over both transmitted zeros and ones. Because of 
this, we group all patterns that are connected via reflection 
or negation into one class of patterns with identical PBER. 
Each class contains at the least two patterns (p and ncg(p) 
because p ^ neg(p), Vp) and at most 4 different patterns 
(p, neg(p), refi(p), and neg(refl(p)) if they are all different) 
and is represented by a unique class index q G 1,...,Q, 
where Q is the number of classes. A labeling can now 
be represented not only by a binary matrix C or by the 
set of pattern indices W, but also by an ordered vector of 
class indices q = [qi, 92, • ■ • , 9m] G {l,---,<9}™, where 
<Zi < 92 < • • • < (Zm- The reason for introducing this vector q 
is that it allows us to easily identify two binary labelings that 
give the same BER. In other words, if two different labelings 
C and C have vectors q and q', resp., they will give the same 
average BER if q = q'. To clarify these definitions, consider 
the following example (see also Example |2]l. 

Example 1 (Patterns for A-PAM): For 4-PAM there ai'e six 
patterns which are grouped into Q = 3 classes as shown in 
Table J] Two of them are RE and four are ARE, as indicated 



Preprint, March S, 2013. 



5 



TABLE I 

Classes of patterns for 4-PAM with their corresponding class 

INDICES q, representatives p, TYPES, DECIMAL REPRESENTATIONS OF 
THE PATTERNS W, VECTORS a DEFINING THE PEER FOR THE ABD, AND 
THRESHOLDS FOR THE REPRESENTATIVES 



<? 


P 


Type 


w 


a 


Thresholds 


1 


[0,0,1,1] 


ARE 


3 12 


[2, 2, 0] 


A = 


2 


[0,1,1,0] 


RE 


6 9 


[4, 2,-2] 


A = -Pi 


3 


[0,1,0, 1] 


ARE 


5 10 


[6,-4, 2] 


= /32 = 



in the third column of Table I] The first column is the index 
of the class and the fourth column contains the decimal 
representations of the indices of the patterns that belong to that 
class. As an example, one of the patterns of a class is shown 
in the second column. These patterns are called representative 
and are used to analyze the performance of the patterns in 
the class. The indices of the representatives are shown with 
boldface in column 4. The fifth column contains vectors a 
defining the PEER for the ABD in (O. The patterns are 
ordered from best to worst PEER for high SNR, as predicted 
by the vectors a. All 4! — 24 valid labeHngs for 4-PAM give 
three vectors q, i.e., three labelings that give different BER: 
q = [1,2] (BRGC), q = [1,3] (NBC), and q = [2,3] (AGC). 
These labelings are listed in the first part of Table HIl together 
with their class indices q, pattern indices W, and vectors a 
defining the BER for the ABD. The labelings are also ordered 
from best to worst BER for high SNR as predicted by the 
vectors a.. 

The next theorem gives closed form expressions for the 
number of classes for length-il/ patterns. 

Theorem 3: For M-PAM, all the length- A/ bit patterns can 
be grouped into Q classes, where 

Q = \{Q%) + if//l)+^''^')^ (15) 

among which Qre classes have only RE patterns, Qare only 
ARE patterns, and Qasy only asymmetric patterns, where 

QHE^iO^ (16) 
Qare = 2'^/^-\ (17) 
QASY = J(a?/g-0-2-/^). (18) 

Proof: Any pattern p must contain M /2 zeros and M/2 
ones, hence, the total number of patterns is equal to (^jyj) ■ 
a pattern to be RE, pi = pM-i+i, i.e., the positions of the M/4 
ones in [pi, . . . ,Pm/2] fully describe the pattern, and thus, the 
number of RE patterns is (^^^4) . There are two members in 
every class of RE patterns, p = refl(p) and neg(refl(p)) = 
neg(p), which gives ( fTSI l. 

For a pattern to be ARE, pi — pM-i+i, i-e., the positions of 
the ones in [pi, . . . ,Pm/2] fully describe the pattern where the 
number of ones in [pi , . . . , pm ^2] is between and M/2. From 
that, it follows that there are 2*^/^ ARE patterns. There are two 
members in every class of ARE patterns {p = neg(refi(p)) 
and refl(p) = neg(p)), which gives (fTTT i. 

All the remaining classes include only asymmetric patterns. 
The number of asymmetric patterns can be obtained by sub- 



TABLE II 

Some common labelings for 4-PAM and 8-PAM with their 

CORRESPONDING CLASS INDICES q, PATTERNS INDICES W, AND VECTORS 
a DEFINING THE BER FOR THE ABD 



M 


Labeling 


q 


W 


a 


4 


BRGC 


[1,2] 


{3,6} 


[6, 4, -2] 


4 


NBC 


[1,3] 


{3,5} 


[8, -2, 2] 


4 


AGC 


[2,3] 


{5,6} 


[10,-2,0] 


8 


BRGC 


[1,2,6] 


{15,60, 102} 


[14,12,-2,0,2,0,-2] 


8 


FBC 


[1,2,10] 


{15,60,90} 


[18,0,4, 10,-8,2,-2] 


8 


NBC 


[1,5,11] 


{15,51,85} 


[22,-4,8,-10,8,-2,2] 


8 


BSGC 


[2,6,9] 


{105, 60, 102} 


[22,10,-8,4,-2,-2,0] 


8 


AGC 


[9, 10, 11] 


{90, 105, 85} 


[36,-18,6,4,-4,-2,2] 



tracting 2(5re and 2(5 are from the total number of patterns, 
i.e., ~ 2(3re — 2(3 ARE- There are four patterns in each 

class, as p ^ refi(p) and p 7^ refl(neg(p)) (or equivalently, 
refl(p) 7^ neg(p)). Using this, ( fTSl ) is obtained. 

Finally, the total number of classes in (flST l is obtained as 

Qre + Qare + Qasy- □ 
Theorem |3] gives the exact number of classes Q for A/-PAM 
constellations. A loose bound on this number of classes was 
previously presented in ||25] eq. (3.14)]. For 4-PAM Theorem|3] 
gives (3 = 3, (5r,e = 1, Qare = 2, and Qasy = 0, which is 
in agreement with Example [T] and Table I] The following two 
examples show the number of classes for 8-PAM and 16-PAM. 

Example 2 (Patterns for 8-PAM): For 8-PAM (M = 8), 
Theorem |3] states that there are Q — 23 classes of patterns, 
Qre = 3 classes of RE patterns and (Jare — 8 classes of 
ARE patterns, each containing two patterns. These classes are 
shown in the first 1 1 rows of Table |lll] and ordered from the 
best to worst PBER for high SNR, as predicted by the vectors 
a. The Qasy — 12 classes of asymmetric patterns, each with 
four members, are listed in the last 12 rows of Table Ulljand 
also ordered in a similar way. By enumerating all 8! = 40320 
valid labelings for 8-PAM, we found 460 vectors q, i.e., for 
8-PAM there exist only 460 labelings that give different BER. 
The classes associated to five of the most well-known labelings 
ai-e: q = [1, 2, 6] (BRGC), q = [1, 2, 10] (FBC), q = [1, 5, 11] 
(NBC), q = [2,6,9] (BSGC), and q = [9,10,11] (AGC). 
These labelings are Usted from best to worst BER in the second 
part of Table HH 

Example 3 (Patterns for 16-PAM): For 16-PAM (M = 16) 
Theorem |3] states that there are Q ~ 3299 classes, (5re = 
35 of the classes contain RE patterns, Qare = 128 classes 
contain ARE patterns, and Qasy = 3136 classes contain ASY 
patterns. The patterns for 16-PAM are not listed in this paper 

V. Thresholds for the BD 

In this section, we show that the thresholds for the BD can 
be found by solving an (A/ — l)th power polynomial equation 
and give a closed-form solutions for 4-PAM and 8-PAM with 
RE and ARE patterns. 

A. Threshold Computation 

The problem of finding the thresholds (3k for the BD 
is equivalent to finding the solutions of l{y) — 0. In the 
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TABLE III 

Classes of patterns for 8-PAM with their corresponding class indices q, representatives p, types, decimal representations of the 

PATTERNS W, VECTORS a DEFINING THE PEER FOR THE ABD, AND THRESHOLDS FOR THE REPRESENTATIVES 



q 


P 


Type 


w 


a 


Thresholds 


1 


[0,0,0,0,1,1,1,1] 


ARE 


15 240 


[ 2, 2, 2, 2, 0, 0, 0] 


ft = 


2 


[0,0,1,1,1,1,0,0] 


RE 


60 195 


[ 4, 4, 2, 2,-2,-2, 0] 


ft = -ft = f{t2) 


3 


[1,1,1,0,1,0,0,0] 


ARE 


23 232 


[ 6,-2, 2, 0, 2, 0, 0] 


ft = -ft = f{ti), ft = 


4 


[0,1,1,1,0,0,0,1] 


ARE 


113 142 


[ 6, 4, 4,-4,-2,-2, 2] 


ft = -ft = f{t2), ft = 


5 


[0,0,1,1,0,0,1,1] 


ARE 


51 204 


[ 6, 6,-4,-4, 2, 2, 01 


ft = -ft = f(t2), ft = 


6 


[0,1,1,0,0,1,1,0] 


RE 


102 153 


[ 8, 6,-6,-4, 4, 2,-2] 


ft = -ft = f{t2), ft = -ft = fit;) 


7 


[0,0, 1,0, 1,0, 1, 1] 

-1-5 v/, -1-5 vy, -L, -Lj 


ARE 


43 212 


[10, -6, 4, -2, 0, 2, 0] 


ft = -ft = /(ts), ft = 


8 


[0,1,0,0,1,1,0,1] 


ARE 


77 178 


[10, 0,-6, 2, 4,-4, 2] 


ft = -ft = f(t2), 

ft = -ft = /(t3), ft = 


9 


[0,1,1,0,1,0,0,1] 


ARE 


105 150 


[10, 0,-4, 6,-4,-2, 2] 


ft = -ft = /(t2), 

ft = -ft = fiu), ft = 


10 


[1,0,1,0,0,1,0,1] 


RE 


90 165 


[12,-6, 0, 6,-6, 4,-2] 


ft = -ft = /(tl), 

ft = -ft = /(ts), ft = -ft = /(t2) 


11 


[0,1,0,1,0,1,0,1] 


ARE 


85 170 


[14,-12,10,-8,6,-4, 2] 


ft = -ft = /(fe), ft = -ft = /(ta), 
ft = -ft = ,f(ti), ft = 


12 


[0,0,0,1,1,1,1,0] 


ASY 


30 120 135 225 


[ 4, 3, 3, 2,-2,-1,-1] 


ft, ft 


13 


[0,0,0,1,1,1,0,1] 


ASY 


29 71 184 226 


[ 6, 1, 2,-3, 1, 0, 1] 


ft, ft, ft 


14 


[0,0,0,1,1,0,1,1] 


ASY 


27 39 216 228 


[ 6, 2,-3, 1, 1, 1, 0] 


ft, ft, ft 


15 


[0,0,1,1,1,0,0,1] 


ASY 


57 99 156 198 


[ 6, 5, 0,-3,-3, 2, 1] 


ft, ft, ft 


16 


[0,0,1,0,1,1,1,0] 


ASY 


46 116 139 209 


[ 8,-1,2,-1, 3,-2,-1] 


ft, ft, ft, ft 


17 


[0,0,1,1,1,0,1,0] 


ASY 


58 92 163 197 


[ 8,-1,3,-2, 2,-1,-1] 


ft, ft, ft, ft 


18 


[0,1,0,0,1,1,1,0] 


ASY 


78 114 141 177 


[ 8, 2,-1,-1,-1,3,-2] 


ft, ft, ft, ft 


19 


[0,0,1,1,0,1,1,0] 


ASY 


54 108 147 201 


[ 8, 3,-6, 3, 3,-2,-1] 


ft, ft, ft, ft 


20 


[0,0,1,0,1,1,0,1] 


ASY 


45 75 180 210 


[10,-3,-3, 6,-4, 1, 1] 


ft, ft, ft, ft, ft 


21 


[0,0,1,1,0,1,0,1] 


ASY 


53 83 172 202 


[10,-3, 1, 0,-2, 1, 1] 


ft, ft, ft, ft, ft 


22 


[0,1,0,1,1,0,0,1] 


ASY 


89 101 154 166 


[10, 0,-3, 1, 1,-3, 2] 


ft, ft, ft, ft, ft 


23 


[0,1,0,1,0,1,1,0] 


ASY 


86 106 149 169 


[12,-6,3,-1,-1, 3,-2] 


ft, ft, ft, ft, ft, ft 



following theorem we show how this can be done for Af-PAM 
constellations. 



Theorem 4: The thresholds for the BD and Af -PAM 
stellations with a pattern p are 



Pk = -^logz„, 
47rf 

where z„ are the real and positive solutions of 

(Af/2-i)(Af/2 + l-i) 

Z = 0, 



and 



M 



-i{y+d(M~2i+l)Y 



(19) 



(20) 



(21) 



Proof: Using Q, l{y) ~ is equivalent to e'*^^-* 
which can be restated as 



(22) 



where the definition of of the A/-PAM symbols was used. 
Factorizing (|22] | gives 



— „-7y^-27dy(M-l)-7d= 



M 



^^g47dy(»-l)g-87d (23) 



Using (EB in (|23]l together with substitution z = e^'t'^v, ^ 
is obtained by setting h{y) ~ and removing the nonzero 
factor preceding the summation in (l23T l. The expression in 
(I20I 1 is an (Af — l)-power polynomial, and thus, it always 
has Af — 1 solutions. Because of the substitution z = e'^'^'^^, 
only the positive (and real) roots need to be considered. □ 

Theorem |4] gives a general expression for finding the 
thresholds for Af-PAM with any pattern p. After finding 
the roots of (|20| |. the thresholds /3k may easily be ob- 
tained from ( fTgj l. The main problem is that finding the 
roots of ( I20I 1 does not admit simple closed-form solutions. 
However, it can always be solved numerically. Fig. |2] illus- 
trates the thresholds for 16-PAM with the pattern P45745 

'Interestingly, for i = 1, . . . , M - 1 with i ^ M/2 and i ^ M/2 + 1, 
the powers of A in )20t are the so-called "triangular numbers". 
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I .516 
I .515 
I .514 
I .513 
I .512 
I .511 
I .510 
I .59 
I .58 



I .54 
I .53 



-10 



10 



20 



7 [dB] 



Fig. 2. Thresholds /3fc in (19) for 16-PAM with the ASY pattern P45745 = 
[1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1] obtained by solving (20} numerically. 
The constellation points are shown with squares. Gray and white areas indicate 
yo and yi, resp. 



[1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1] obtained by solving ^ 
numerically. In the following two sections we show how this 
problem can be solved analytically for 4-PAM and 8-PAM 
with RE or ARE patterns. 



B. Thresholds for 4-PAM 

The following theorem shows how the thresholds are found 
for 4-PAM with any pattern. 

Theorem 5: The thresholds /3fc for any pattern for 4-PAM, 
as listed in the last column in Table U can be expressed as 



/3i = -/33, 



1+P1P4^+V/(1+P1P4V4)2 - 4A2 



2A 



(24) 
(25) 

(26) 



where A is given by (ISTT i. 

Proof: The proof is given in Appendix |B] □ 
Theorem |5] gives closed-form expressions for any pattern 
for 4-PAM, and thus, it allows us to compute the BER for 
all the labelings in the first part of Table HI] The results in 
Theorem |5] can be shown to coincide to those in |4, eq. (10)] 
when the BRGC is considered. 



1, 2, . . . , 11 for 8-PAM can be expressed as 

= —P&-k 

1 




with 

t2 = 



T + + 



6piA3 
1 

6piA3 
1 



B 



A is given by (1211 1. 

T ^2{psA' - P2), 



^(1 + V3j)C _ l-VSj 

B ^ 
^(1 - Vi])C 1 + V3j 
B ^ 



B 



B 



(27) 
(28) 

(29) 
(30) 
(31) 



B^\J^D-^-AC^-PiPsD, 

C^7A6+p2P8A3-3piP3A + l, 

D ^7piA9 - 12pip2P8A6 - 18p3^* 

+ 3]5i(l + 9p4P8)^^ - 9p2P3P8^ + 2piP2P8, 

and the relationship between k and n in dZTl l for the different 
classes q is listed in the last column of Table [III] As the 
relationship between n and k depends on the particular pattern, 
the representative of the class (the second column of Table HiH) 
should be used in the presented equations. 

Proof: The proof is given in Appendix |C] □ 
Theorem |6] shows how to analytically obtain the thresholds 
for the BD with 8-PAM and any RE or ARE pattern, for 
instance, thresholds shown in Fig. [T] Using these results, 
the PBER can be calculated using Q, which gives PBER 
expressions for 11 out of 23 classes, or equivalently, for 56 
different labelings, including the 5 shown in the second part 
of Table In] 

Remark 2: In the high SNR regime, i.e., 7 — 00, all the 
thresholds in Theorems |5] and |6] tend to midpoints, i.e., the 
same constant thresholds used in the ABD for all SNR. This 
fact can easily be proven analytically for 4-PAM by evaluating 
lim^_>.oo and applying I'Hopital's rule. For 8-PAM a similar 
proof exists, however, in this case it is not straightforward due 
to the complexity of the threshold expressions. These results 
can be intuitively understood from the fact that the max-log 
approximation in (JSjl becomes more precise when the SNR 
increases, and hence, the thresholds for the BD and ABD are 
expected to coincide when 7 cx). 



C. Thresholds for 8-PAM 

The following theorem shows how the thresholds may be 
found for 8-PAM with RE or ARE patterns. These patterns are 
of great value because all the most commonly studied labelings 
(e.g., BRGC, NBC, FBC, BSGC, AGC) can be composed from 
them (cf. Table HB. 

Theorem 6: The thresholds f3k for the patterns in the classes 



VI. Numerical Results 

In Fig. [3] we show the thresholds given by Theorem |6] for 
the pattern pj^gg {q = 10) for 8-PAM. The figure is symmetric 
with respect to zero due to the symmetry of the pattern. At 
7 « 5.3 dB the pairs of thresholds (32 and /33, and f3c^ and /Jg 
merge and become virtual for all 7 < 5.3 dB. All the virtual 
thresholds shown with dashed lines satisfy the conditions in 
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7 [dB] 



Fig. 3. The thresholds for p^gg = [1, 0, 1, 0, 0, 1, 0, 1] (q = 10) for 8-PAM 
in Theorem |6] vs. SNR. Virtual thresholds are shown with dashed lines. The 
thresholds for the ABD fi^ and the constellation points (squares) are also 
shown. Gray and white areas indicate and 3^1. resp. 



Theorem |2] As expected (see Remark |2), when 7 — > cx), the 
BD thresholds coincide with the ABD thresholds. 

The PBER for 8-PAM with some selected patterns from Ta- 
ble ini] using (|9]) is presented in Fig. |4] The thresholds are 
calculated analytically for q = 3, 10 and numerically for 
q = 16, 22. For very low SNR the gap between the BD and the 
ABD can reach up to several dB, however, this gap decreases 
when the SNR increases. The same conclusion can be drawn 
for all other patterns except for q = 1, as in this case only one 
threshold exists /3j\//2 — for all SNR. Hence, for q = 1 the 
BD and the ABD have the same performance for A/-PAM. To 
conclude, we present in Fig. |5] the BER for 8-PAM with the 
labelings in Table [III From the presented results we conclude 
that the BD outperforms the ABD, however, for any BER of 
practical interest (below 0.1), the difference between the BD 
and the ABD is negligible. 

VII. Conclusions 

We proposed a general approach for estimating the per- 
formance of the optimal bit-wise demodulator and presented 
closed-form expressions for the BER for 4-PAM and 8- 
PAM with different labelings. We conclude that a suboptimal 
symbol-wise demodulator shows no loss compared to the 
optimal demodulator for all the SNR of interest, which justifies 
its use in practical systems. 

The derived BER expressions for the optimal demodulator 
can be used to calculate the mutual information (MI) of BICM 
when the demodulator makes hard decisions on the bits. The 
optimal bit-wise demodulator does not necessarily maximize 
the MI. Finding the hard-decision demodulator that maximizes 
the MI is left for future work. 

The proposed technique for finding the zero crossings of the 
L-values for 8-PAM works only for reflected or anti-reflected 



0.5 



oa 

Oh 

10-' 





j~" "~ □ - ^ 
















<j = 3 

* g = 10 
□ q = 16 
V q = 22 













5 10 15 

7 [dB] 



Fig. 4. The PBER for 8-PAM with ARE (<j = 3) and RE (q = 10) patterns 
and ASY patterns (q = 16, 22). Solid lines correspond to the BD and dashed 
lines to the ABD. The threshold for the BD were obtained using Theorem |6] 
for g = 3, 10 (the thresholds for g = 10 are shown in Fig. [5] and solving )20t 
numerically for q = 16, 22. 



0.5 




7 [dB] 

Fig. 5. The BER for 8-PAM with the 5 labelings in Table Solid lines 
correspond to the BD and dashed Unes to the ABD. 



patterns, which includes 11 out of 23 classes of patterns. 
Extending these results to the remaining classes of patterns for 
8-PAM is left for further investigation as well as generalizing 
the results to arbitrary M. 
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Appendix A 
Proof of Theorem[T] 

Let Vi,k be the conditional probabilities 

v^,i=Pr{Y <l3i\X = s,} 

= 1-Q((/3i-5,)v/2^), (32) 
v,,k = Pr{f3k-i<Y <l3k\X ^s,} 

= Q ((/3fc_i - s,)y2^) - Q - s,)V^) , (33) 

V^,M = Mhl^l < Y\X = S,} 

= Q((/3Af-i-sOV27) , (34) 

where /3fc for fc G /C are the thresholds and none of them is 
virtual. The PEER in ([8]) can now be rewritten as 



Solving h{z) = gives the three roots zi ~ —pipi and 



^ M M 



(35) 



i=l k=l 



where e^^fc =pi®pj,e {0, 1}. 

Using (I32]l-(l34ll the PEER in ((351) can be expressed as 



P 



M 



M M 



■1=1 1=1 fe=2 



i=l k=l 
M A/-1 



i=l fc=l 



(36) 



where = = M/2 was used. To obtain 

the expression in (|9]l, we express ei,k+i — Ci^k in (l36] l as 

ei,fc+i - ei,fc = pfe+i ® Pi - _pfc ® Pi (37) 

= (pfe+i -Pfe)(l-2p0, (38) 

where the identity pi = pipj + piPj was used together 
with pi = l -pi. 



Appendix E 
Proof of Theorem[5] 

Define the function h{z) as 

h{z) = p4Az^ + p3z'^ + p2Z + piA. (39) 

According to Theorem |4] for 4-PAM with a pattern p = 
[pi,P2,P3,P4], equation h{z) ~ needs to be solved in order 
to find the thresholds. The patterns for 4-PAM are either RE 

(Pt = pM+i-i, Vi) or ARE {pi = ~Pm+i-i, Vi). Therefore 



h{z) = piAz^ + p2Z^ ± p2Z ± piA, 



(40) 



where the upper and the lower signs correspond to RE and 
ARE patterns, resp. Using p| = 1 and the fact that pipi = ±1 
for RE and ARE patterns, resp., h{z) can be factorized as 

h{z) = Pi{z + pipi){Az'^ + (pij52 - APiP4)z + A). (41) 



PiPaA - piP2 ± V {PiPaA - PlP2)^ - 



2A 



(42) 



For q — 1 (where pipA = — 1 and piP2 ~ 1) the root 
zi — 1, that used in ( fT9l l gives the threshold (32 = 0. The 
other two roots in (|42] | are complex for all SNR and do not 
result in thresholds. 

In a similar way, for q ~ 3 (where pij52 — —1 and 
pipA = —I) (32 =0. When A < 1/3, or equivalently, when 
7 > 51og3/8 « —1.63 dE, the roots in (l42l i are real and 
positive resulting in thresholds /3i and /S^ by using ( fT9l l. When 
A > 1/3 (low SNR), the roots in (l42l l are complex and can 
no longer be used in ( fT9] l for calculating the thresholds. To 
overcome this, \z2\ and jzaj are used in the calculation of the 
thresholds, which together with piP2 = ~1 gives ( |26] |. The use 
of I • I does not affect the result when the roots are real. When 
the roots are complex, their absolute values are equal to one 
(can be seen from (|42]|). which gives two virtual thresholds 
Pi ~ (3^ ~ merging with the zero-threshold 132 at around 

— 1.63 dE. Theorem 12] allows the use of these thresholds in 
the calculation of the PEER. 

Finally, for q — 2 (where piP2 = — 1 and j5ip4 = 1) zi = 

— 1, which results in no threshold. The two roots in (l42t are 
positive for all SNR, resulting in the thresholds /3i and /Ss by 
using iT% . As the roots are positive the use of | • | does not 
affect the result, which gives ( |26] |. This completes the proof. 



Appendix C 
Proof of Theorem[6] 

Define the function h(z) as 

h{z) = psA^z'^ + PtA^z^ + peAz^ + p^z* + p^z^ 

+ P^Az^ + P2A^ z + PiA^ . (43) 

According to Theorem |4] for 8-PAM with pattern p ~ 
[pi,P2, ■ ■ ■ iPs] equation h{z) = needs to be solved in 
order to find the thresholds. For RE and ARE patterns, 
pi = zLpM+i-i, Vi, where the upper and the lower signs 
correspond to RE and ARE patterns, resp. Using this property, 
h{z) for RE and ARE patterns is 

h{z) = PiA'^z^ + P2A'z'' + PsAz^ + Paz^ 

±p4z^ ±p3Az'^ ±p2A^z±piA'^. (44) 

Factorizing (l44l i gives 

h{z) = (z ± 1) (piA^e + [P2A''tPiA^] 

+ [PiA^ T P2A^ + psA] z* + [tPiA^ + p2A^ T P3A+ p^] z^ 

+ [piA"^ T P2A' +P3A] z^ + [P2A^ T Pi A"] z + . 

(45) 

Rearranging the terms in (l45t h{z) can be written as 

h{z) = z^{z ± 1) {piA^iz^ + z-^) + [P2A^ T PiA^] (2^ + z-^) 
+ [piA^ T p2A^ + P3A] {z^ + z-^) 

+ [tPiA^ + p2A^ T P3A + P4] ) • (46) 
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Using the substitution 



t{z) = z + z \ 



(47) 



( |46] l can be modified to 

hiz) = z^{z±l) {piA'^tizf + [p2A^TPiA'^] t{zf 
+ [-2PiA'^TP2A^ t{z) 

+ [±PiA'' -p2A'TPzA+Pi\). (48) 

Finding positive roots of h{z) — can now be done 
analytically. For ARE patterns the second factor in (@8) gives 
a root equal to one resulting in /34 — 0. The roots of the last 
factor in (l48T l need to be found. As a first step we solve it 
with respect to t{z), where the roots t„ are shown in (|29]l- 
(|3TT l. where ±1 was replaced by +pips to distinguish between 
RE and ARE patterns. As z should be real and positive, only 
real and positive t„ need to be considered. Two out of three 
roots tn may combine into a complex conjugated couple, but 
the third root is always real. Every positive root tn gives two 
roots for z in (|45T l, which can be found from ( |47] i as 



^2n-1.2n 



(49) 



where Z2n-i = l/z2n- Due to ( fT9l ) and the symmetry of the 
patterns, these two roots give the two thresholds /3fe = — /Jg-fe, 
which justifies the first equality in ( |27] |. When t„ is real and 
tn > 2, the roots in (|49] | are positive and give thresholds 
(3k — — /?8-fc- Because of i„ is real and the roots Z2„-i and 
Z2n are positive, the use of | • | (three times) in (l28l l does not 
change the result. By analyzing all the roots i„, the thresholds 
were found and listed in Table HID where the last column shows 
the relation between the threshold /3k and the roots t„ shown 
in dal-dlB- 

For the listed thresholds in Table HUl t„ in (|29ll-(|3B is 
never real and negative, however t„ can be either complex or 
real with < t„ < 2 for some 7 < 70, resulting in virtual 
thresholds /3k and /^s-fc- In what follows, we show that these 
thresholds are equal to each other when using (|28^ . i.e., they 
fulfill the conditions in Theorem |2] First, consider the case 
when tn is real but < i„ < 2. In this case 2:2n-i and 
Z2n are complex with unit magnitude and according to (l27T l. 
the corresponding thresholds f3k — —/Ss-fe are equal to zero. 
By analyzing the thresholds for all the RE and ARE patterns, 
we find that /Ss = —/S^ for q — 3,6,9,11 and /3q = —P2 
for (7 = 5, which are separated by either no threshold or by 
the threshold /34 = 0. These thresholds can be used in (|9]l 
according to Theorem |2] 

Second, consider the case when tn is complex. One of the 
other two roots of (|48] | tn' giving f3k' = —/3s-k' is such that 
tn' = tn*, which means that |i„| — This leads to two 
pairs of the thresholds f3k — /3k' and /^s-fc = Ps,-k'- Revising 
the thresholds for all the RE and ARE patterns we conclude 
that corresponding thresholds are: /3q = —/32 and /Jy = — /3i 
for g = 8, 11 and ^5 = -/Jg and /3g ^ -^2 for q 7, 10. 
These thresholds can be used in (|9]l according to Theorem |2] 
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